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We study theoretically the nonstationary circuit QED system in which the artificial atom tran- 
sition frequency, or the atom-cavity coupling, have a small periodic time modulation, prescribed 
externally. The system formed by the atom coupled to a single cavity mode is described by the 
Rabi Hamiltonian. We show that, in the dispersive regime, when the modulation periodicity is 
tuned to the 'resonances', the system dynamics presents the dynamical Casimir effect, resonant 
Jaynes-Cummings or resonant Anti-Jaynes-Cummings behaviors, and it can be described by the 
corresponding effective Hamiltonians. In the resonant atom-cavity regime and under the resonant 
modulation, the dynamics is similar to the one occurring for a stationary two-level atom in a vi- 
brating cavity, and an entangled state with two photons can be created from vacuum. Moreover, 
we consider the situation in which the atom-cavity coupling, the atomic frequency, or both have 
a small nonperiodic time modulation, and show that photons can be created from vacuum in the 
dispersive regime. Therefore, an analog of the dynamical Casimir effect can be simulated in circuit 
QED, and several photons, as well as entangled states, can be generated from vacuum due to the 
anti-rotating term in the Rabi Hamiltonian. 
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I. INTRODUCTION 



Over the last few decades nonstationary processes in cavities have received considerable attention. One such process 
is the nonstationary or dynamical Casimir Effect (DCE) - in particular, the creation of photons from vacuum, or 
another initial field state, in cavity whose geometry [H, or material properties (1, IE, S 0, B [§| have a periodic 

time dependence, with the modulation frequency close to twice the unperturbed field eigenfrequency. Nowadays, 
DCE in cavities is a well studied problem, with a variety of theoretical predictions concerning the number and the 
statistics of created photons, as well as the influence of detuning, dissipation, boundary conditions, geometry and 
nonperiodicity of the modulation (see [TO, [Hi [l^l for an extensive list of references) . To date, DCE has not been 
observed in laboratory, however several concrete proposals have app eared over the last years [H, [H, [H, [H, [13, [3 , 
with some of them being currently implemented experimentally |l9l . [20l . [2]| . 

The interest in nonstationary processes in cavities rea ppe ared over the last 5 years due to the progress in the field 
of Cavity Quantum Electrodynamics (cavity QED p^ . |23| | ) in the condensed matter systems, e.g. semiconductor 
quantum dots [H, [H, [H, [23, [H, [2^ , polar molecules [30l |3l| and superconducting circuits [H, [13, [11] coupled to 
high-Q resonators, the latest architecture known as circuit QED [H, [35'|. In cavity QED, the effective 2-level atom 
is coupled to the field inside the resonator via the dipolc interaction, allowing for observation of the light-matter 
interaction at the level of single photons and single atoms. The new ingredient in the solid state cavity QED, in 
particular circuit QED, is the possibility of engineering and manipulating the properties of the artificial 2-level atom 
and the resonator 33, 36, 31^, as well as the interaction strength between them [38], either during the fabrication, or 
in situ. 

Recently, the strong resonant and the strong dispersive coupling limits between the artificial atom and a single 



cavity mode were observed experimentally in circuit QED [33, |39| and other solid state architectures [2J, [25j, [261, 
[27l . [29j . Moreover, the single photon source [13], single artificial-atom maser [4lj, multiphoton Fock states 42] 
and interaction between two artificial atoms (qubits) (stI . [43| were implemented experimentally, among many other 
important achievements (see [llj for more references) . Besides, the circuit QED architecture benefits from robust 
readout schemes of the atomic and the resonator states [ss", [3, [43, [4^ , relatively low dissipative losses 

[sot , state preparation techniques [49i] and real time manipulation of the atomic transition frequency via electric and 
magnetic fields [H, [13, [12] or nonresonant microwave fields [H, \^ . 

Harnessing the tunability of the atomic transition frequency, the authors in [sol [5l| proposed an experimental 
implementation of the Landau-Zener sweeps in circuit QED, when the atomic frequency increases linearly in time. This 
allows for generation of single photons, as well as cavity-atom entangled states. In a similar direction, the generation of 
the quantum vacuum radiation by modulating periodically the vacuum Rabi frequency of an intersubband transition 
in a doped quantum well system embedded in a planar microcavity was considered in [TtI . [T^ [s^ . [s^ , and emission 
of photon pairs was predicted for 'resonant' modulation frequencies. A preliminary theoretical study of the feasibility 
of realizing the dynamical Casimir effect with a quantum flux qubit in superconducting quantum nanocircuits, as 
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well as the detection of the generated photons, was reported in [5j|. On the other hand, the freedom of controlling 
in real-time the atomic frequency is currently being explored to couple/decouple one or several qubits to/from the 
cavity mode in order to implement quantum lo gic o perations [stI. [43L l48l| . 

In this paper, following the original proposal 55j, we study nonstationary processes in the solid state cavity QED, 
where a single cavity mode is coupled to a single artificial atom whose transition frequency, or atom-cavity coupling, 
have a small periodic time modulation, prescribed externally. Such a control over the transition frequency, with 
compatible modulation frequency, can be achieved in circuit QED with present or near- future technology [stI [43 . |48| . 
We show that in the dispersive atom-cavity regime and under the resonance conditions one obtains completely different 
effective regimes of the atom-field interaction, which can be approximately described by the resonant Anti-Jaynes- 
Cummings (AJC), resonant Jaynes-Cummings (JC) or the dynamical Casimir effect (DCE) Hamiltonians. Moreover, 
in the resonant atom-cavity regime, the system dynamics resembles the behavior of the DCE in the vibrating cavity 
containing a resonant two-level atom, and entangled states with at most two photons can be generated from vacuum. 

We also consider the case, in which the atomic frequency, the atom-cavity coupling parameter, or both have a 
small nonperiodic time modulation, prescribed externally. Namely, we suppose that the modulation is given by the 
sum of two harmonic functions with different amplitudes and frequencies. We deduce an effective Hamiltonian in the 
dispersive regime, and show that photon generation from vacuum is possible for fine tuned modulation frequencies. 

Thus, we demonstrate the possibility of simulating the DCE in circuit QED using a single nonstationary atom, 
instead of a macroscopic dielectric medium as in Refs. 0, B 0- As applications, it could be possible to create 
excitations, either photonic or atomic, from the initial vacuum state \g,0), generate nonclassical states of light and 
realize transitions between the states {\g,n), \e,n± 1)} in the dispersive regime, where \g) and |e) are the atomic 
ground and excited states, respectively, and \n) is the Fock state of the cavity field. A related problem was recently 
studied in (56j . where it was suggested that lasing behavior and the creation of a highly non-thermal population of 
the oscillator, as well as the cooling, could be implemented using an analogous scheme in the near future. 



II. NONSTATIONARY CIRCUIT QED WITH PERIODIC MODULATION 



We assume that the atomic transition frequency Q (t) is given by 

n (t) = no + eft . (1) 

Here flo denotes the bare atomic frequency, e is a small modulation amplitude, e ^ Oq, and ft is an arbitrary periodic 
function of time, prescribed externally 

oo 

/t = ^ [sfc sin fcyyt -|- Ck cos krjt] , (2) 

k=0 

where rj is the modulation frequency and {sk,Ck} form a set of coefficients describing the time modulation. The 
cavity frequency uj and the atom-cavity coupling parameter go are constant, so at the charge degeneracy point [38| 
the system is described by the Rabi Hamiltonian (RH) 

H = Ho{t)+gn{a + a^){a++<j^), (3) 

with a (at) being the cavity annihilation (creation) operator, cr+ = \^){g\ s-nd ct_ — |5)(e|. The free Hamiltonian is 

Ho{t)^con+^a,, (4) 

where n = a^a is the photon number operator, cr^ ~ \^){^\ ^ \9){9\ ^-nd. we assume h — 1. In the stationary case, 
e = 0, one can perform the Rotating Wave Approximation (RWA) and obtain the standard Jaynes-Cummings (JC) 
Hamiltonian [Hjl, which has been verified in several experiments over the last few years [H, [H, [H, [iOl . However, 
in the nonstationary case, as well as under strong dephasing noise [Qlj, the anti-rotating term (atj- -I- a^cr^) cannot 
always be eliminated. Moreover, it is responsible for producing an analog of the DCE and creating photonic and 
atomic excitations from vacuum under modulation 'resonance' conditions, as shown below. 
In the interaction picture with respect to Hq (t) the interaction Hamiltonian reads 

Hj = go (e*""ao'+ + e*"+a^o'+ + h.c.) , (5) 

where h.c. stands for the Hermitian conjugate and 



/ dT[n{T)±i 

Jo 



(6) 
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All the information about the influence of the external modulation on the system dynamics is contained in the time- 
dependent coefficients exp(zS±), which can be significantly simplified by tuning the modulation frequency 77 to the 
'resonances'. We have explicitly 



00 



fe=i 



where we defined the complex coupling constant 



and parameters 



g = go exp 



At 



Sk 



k=l 



Ck + isk 



2k 



(7) 



(8) 



(9) 



(10) 



III. ANTI-JAYNES-CUMMINGS (AJC) RESONANCE 

The 'Anti-Jaynes-Cummings' (AJC) resonance occurs for 

V = Vajc = A+ - ^, 



(11) 



where |^| ^ A_). is a small 'resonance shift' to be adjusted afterwards. Assuming a reasonable experimental condition 
e/?7 <C 1, we expand g^ exp(iS-i-) in ([5]) to the first order in e/rj and make RWA, obtaining 



where the modulation induced dimensionless coupling is 



e = Ai- 



(12) 



(13) 



A. Dispersive regime 

In the dispersive regime, go\/{n}/ |^-| ^ 1; where (n) is the mean photon number, the Hamiltonian (|12p can be 
approximated by [s^l 



H 



(1) 



where 



and 



9e*«*a V+ + h.c.) +S{n + 1/2) a. 



6 = So+0{g'o/A+) 



(14) 
(15) 

(16) 



is the standard dispersive shift 57]. After the unitary time-dependent transformation defined by the operator 

U = exp [i^ta^/2] (17) 

we get the effective AJC Hamiltonian 



Hajc = - ^ i±i(l±2n) ^ ^ 



(18) 
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For the field initially in the Fock state, by adjusting the resonance shift ^ in order to make |^ + 5 (1 + 2h)\ small 
compared to \gO\, one obtains approximately the resonant AJC Hamiltonian 

H^^j(j ~ gea)a+ + g*e*aa^. (19) 

From the physical point of view, the external modulation supplies the energy w + flp necessary to create one photon 
and one atomic excitation simultaneously. Thus, one can create the superposition of states |e,l) and |.g,0) starting 
from the initial vacuum state |(7,0) and fine tuning the resonance shift. 

In figure [T] we show the exact dynamics for the AJC resonance, obtained through the numerical integration of the 
Schrodinger equation with the initial Hamiltonian Q, setting the experimental circuit QED parameters fJo/ti' = 1.4, 
gQ/oj = 2 • 10^^ and assuming s/lo — 0.2. In all the simulations throughout this paper we consider the harmonic 
modulation, ft = sin rjt. In the figure [1^ we plot the number of created photons 

N={n)^ no, (20) 

for the initial Fock state \g, tiq) as function of the dimensionless time |^| g^t , assuming the resonance shift ^ = 
—26o {1 + no). In figure [T|d we plot the atom excitation probability Pp. From figures [Iji and[T|D one can observe the 
simultaneous generation of one photon and one atomic excitation, in agreement with the effective AJC Hamiltonian 
(lisp . However, the maximal mean number of created photons is slightly less than 1 and depends on the initial photon 
number no, meaning that the resonance shift was not adjusted precisely. In figure [TJ; we plot for the initial number 
state \g, no — 5) for different resonance shifts ^ = —2So {1 + no + x), x = 0,1, 2, 3, 4, demonstrating that it is possible 
to find an optimum ^ which reproduces the resonant AJC Hamiltonian and allows for the creation of one photon and 
the atomic excitation simultaneously. 




FIG. 1: The results of the numerical integration of the Schrodinger equation with the Rabi Hamiltonian ([3]) for the modulation 
frequency rj = A+ + 25o (1 + no). We see an AJC behavior, as given by the effective Hamiltonian (|18p . a) The number of 
created photons N and b) the atom excitation probability Pe versus the dimensionless time for the initial Fock states \g,nQ}, 
with no — (line 1), no = 3 (line 2) and no = 6 (line 3). c) The dependence of N on the small change x in the resonance shift 
^ = —5 (1 + no + x) for the initial number state \g, no = 5). For \0\ got from 2 to 2.3 the curves correspond to a; = 4, 0, 3, 1, 2, 
from below. By adjusting x one can optimize the photon generation process. 



B. Resonant regime 

In the resonant regime, |A_| /go <C 1, the effective Hamiltonian method is not applicable. So we apply the method 
of slowly varying amplitudes on the Hamiltonian (|12p , repeating the procedure [631 employed originally for studying 
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the photon generation from vacuum due to the DCE in a vibrating cavity containing a resonant (stationary) two-level 
atom (63l [m . [g^. We find that for the initial state \g,0) the photon generation occurs for two possible resonance 
shifts 

^± = A_/2±V2.go (21) 

and one gets the following non-zero probabilities Pan, with a ~ (e,.<7) denoting the atomic state and n the photon 
number 

Pgo « cos2(xi) 

Pel - sm^{y + q)sm^xt) (22) 
Pg2 « cos^ {y - q) sin^ (xt) • 

Here 

X^\g0\sm{y + q), tany « (^^^^±|^ j ^ (23) 

and g = (7r/2) for f_ (^+). 




FIG. 2: The time dependence of the probabilities Pgo, Pel and Pg2, obtained by the numerical integration of the Schrodinger 
equation with the Rabi Hamiltonian ^ in the resonant atom-cavity regime for the AJC resonance with properly adjusted 
frequency shift, a) A_ = 0, ^ = ^_ and e/ujr ~ 2- 10~^. b) A_ = W~^go, ^ = C- s-nd e/ajr = 10~^. Not more than 2 photons 
can be generated from vacuum, and the atom excitation probability is limited by the value ~ 1/2, in agreement with equation 
((22)). 

We illustrate this behavior in figure[2l where we show Pgo, Pel and Pg2 in the resonant atom-cavity regime under the 
AJC resonance. In figure[2^ we consider A_ = and ^ = using the parameters ^o/f-^ = 4- 10~^ and e/cu — 2- 10^^. 
In figure[2lD we set A_ — 10~^5o, C — C-i go/^^ ~ 4 • 10^^ and e/uj — 10"^. The results are in excellent agreement 
with the theoretical predictions, equation P^ . demonstrating that a superposition of states |5,0), |e, 1) and |.g,2) is 
created from the initial vacuum state \g, 0). This dynamics is similar to the one studied in the context of DCE, where 
a resonant (stationary) two-level atom is fixed inside a cavity whose boundary is oscillating with the frequency close 
to 2u! [sll- In both cases not more than two photons can be created from the vacuum state \g, 0) and the probability 
of exciting the atom is limited by the value 1/2. This similarity is not surprising, because for the AJC resonance the 
modulation frequency is 77 « A+ « 2a;, and the artificial atom plays the role of the stationary two-level atom and the 
cavity modulating mechanism at the same time. 
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IV. JAYNES-CUMMINGS (JC) RESONANCE 



In the dispersive regime the 'Jaynes-Cummings' (JC) resonance occurs for 
For A„ > 0, repeating the steps that led us to p^ . we get the effective JC Hamihonian 



H 



JC 



-<^z + {gS aa+ + h.c.) . 



(24) 



(25) 



ff A_ < 0, we obtain upon the replacements 9 —9* and ^ Thus, by employing the JC resonance and 

adjusting the resonance shift ^, one can couple the subspaces {\g,n), \e,n— 1)} when the atom and the field are far 
off resonant, since the external modulation supplies the energy difference jfio ~ ^| necessary to couple the atom and 
the cavity field. 
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FIG. 3: The results of the numerical integration of the Schrodinger equation with the Rabi Hamiltonian Q for the JC resonance, 
rj — A- + 25o (1 + no), considering initial number states \g,no} with no — 1 (line 1), no = 3 (line 2) and no — 7 (line 3). 
a) The number of created photons A'^ and b) the probability of finding the atom in the ground state Pg as functions of the 
dimensionless time demonstrate a transfer between the states \g,no) and |e,no — 1) in the dispersive regime, c) for the 
initial number state |(7,4) and three lower JC resonances (?) = t)!^', K = 1,2,3 as shown, see section IVl)) for the resonance 
shift ^ = -lOSo- 

The JC behavior is illustrated in figure [3] for the initial Fock state \g,no). In figures [3^ and [3)3 we plot A'' and 
Pg = 1 — Pe, respectively, using the same parameters as in figure [1] and setting 77 = A_ + 25 (1 + no). As expected 
from the effective JC Hamiltonian one has transitions between the states {\g,n), \e,n— 1)} for large detuning, 
|A_| ^ g. Thus, in the dispersive regime one could engineer entangled states with several photons from the initial 
vacuum state \g,0) by alternating between the AJC and JC resonances and controlling the time interval of each 
resonance. 



V. DYNAMICAL CASIMIR EFFECT (DCE) RESONANCE 

In the dispersive regime, the 'dynamical Casimir effect' (DCE) resonance occurs for 

V = Voce = A+ - A_ - 2^ = 2^^ - 2^. 



(26) 
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Performing RWA in the interaction Hamiltonian ([5]) we obtain to the first order in e/77 

Hi ~ ge'^-'aa+ + g0e'('^-+2«)*aV+ + h.c. (27) 

Applying the unitary time-dependent transformation defined by the operator 

Ui = exp {i [in + (A_ + (Jzl'A ^1 , (28) 

we obtain the time-independent Hamiltonian consisting of the JC Hamiltonian plus the anti-rotating term multiplied 
by the adjustable coupling 176' 

Hf^ ~ifi+ ^ a^ + (ffag+ + g0a^(T+ + h.c.) . (29) 

Since jg/A-j <C 1, we can obtain an effective Hamiltonian by applying a sequence of small unitary transformations 
on Hj^^ [6^ and performing the Hausdorff expansion in some small parameter at each step. Here our expansion 
parameter is e = jg/A-j and, assuming that 9 O (e), we consider terms up to the second order in e. Applying first 
the 'rotating' unitary transformation defined by the time-independent operator 

Ur = exp [{gaa+ - h.c.) /A_] (30) 

followed by the 'antirotating' transformation defined by the operator 

Ua = exp [(.90aV+ - h.c) / (A_ + 2?)] , (31) 

we obtain the effective Hamiltonian 

H.ff^UaUrH^^UlUl (32) 

of the following form (considering terms up to the second order in e) 

Heff ~{i + 5(Tz) h + 5a, {0*a'^ + h.c.) + ^^^^^cr, - ^ igana+ + h.c.) + O {e^) . (33) 



Finally, moving to a rotating frame by means of the unitary time- dependent transformation defined by the operator 

U2 = exp [~i (A_ + <5 + a,t/2] , (34) 

we obtain 

Hdce ^ (e + Sa,) fi + 5a,{ea^^ + e*a^) - — {ge'^-'ana+ + h.c.) + O (e^) , (35) 



where we neglected {5 -\- i) compared to A_ in the exponential of the third term. 

This is the main result of this section. The first two terms of the effective Hamiltonian ([35]) form the DCE part and 
the remaining terms represent corrections, oscillating with a high frequency ~ A_. For small (n) <^ Uc ^ (A_/(7) , 
the contribution of these corrections is relatively small and a, becomes approximately a constant. Assuming that 
initially the atom is in the ground state, then for (fi) <g; ric one has cr^ « — 1 and the right-hand side of equation (|35p 
transforms into the DCE Hamiltonian [j, [T]| 

HDCE^{i-5)n~6{ea)^ +9*0^). (36) 

If the atom is initially in the excited state, a similar effective Hamiltonian is obtained through the substitution 5 — > —5. 
Therefore, by adjusting the frequency shift to ^ = ±J, depending on the initial atomic state [68| . we have photon 
pairs creation from vacuum, as well as field amplification, due to an analog of the DCE. 

In this case the DCE is simulated by the atomic transition frequency modulation through the strong atom-cavity 
coupling. However, the photon generation process is not steady because after several photons have been created the 
third and higher order terms become important and the photon generation is interrupted. Indeed, the third term 
describes the non- resonant atomic excitation by means of the absorption of photons from the cavity. Nevertheless, 
the Hamiltonian (j35|) shows that it is possible to simulate DCE and generate several photons from vacuum using a 
single artificial atom. 
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A. Qualitative explanation and numerical results 



This phenomenon can be quahtatively understood as follows. In the dispersive regime the atom acts as an effective 
nonlinear capacitance [38| . pulling the cavity frequency to approximately 

u;{t)^L, + -^^^^{u; + aJ)~aJ^ft. (37) 
A_ + ejt A_ 

Consequently, it is ex pec ted that the periodic modulation of ft with the modulation frequency close to 77 « 2 (lj ± (5) 
will lead to DCE [Til. l62j| . for which the photons are generated as long as the modulation is present. The energy 2cj 
necessary to create pairs of photons is provided through the atomic frequency modulation and the resonance shift ^ 
must be adjusted [g^I in order to get the constructive interference on the cavity field [HI, [HI. However, in our case 
the atom gets entangled with the field due to the third term in (|35p . by which the atom acquires some probability of 
being excited through photon absorption and the photon generation process cannot continue asymptotically due to 
the loss of constructive interference [the first term in (j35p becomes non-zero]. This is different from the usual DCE 
situation, in which theproperties of the macroscopic linear, lossless and non-dispersive dielectric medium inside the 
cavity are modulated [sHMUl s^-nd the field mode does not get entangled with individual atoms. 




FIG. 4: The initial stage of the time evolution of a) the number of created photons A'' and b) the atom excitation probability 
Pe, obtained by the numerical integration of the Schrodinger equation with the Rabi Hamiltonian Q for rj — 2 {ur — So) , for 
the initial Fock states \g,no) with no = (line 1), no = 1 (line 2) and coherent states \g,a) with \a\^ — 1/2 (line 3) and |a|^ — 1 
(line 4). There is photon creation due to an analog of the dynamical Casimir effect, as described by the first two terms in 
the effective Hamiltonian (|35|l . c) A'' and d) Pe for larger times and initial state \g,Q). Here = 2 {cor — xSo) with x = 0.95 
(line 1), X — 1.25 (line 2), x = 1.15 (line 3), j: = 1 (line 4) and x = 1.05 (line 5). The photon creation is very sensitive to the 
resonance shift, and the generation process is interrupted due to the atom-field entanglement, as expected from the third term 
in equation (|35p . 

In figure S] we present the results of numerical integration of the Schrodinger equation with the Rabi Hamiltonian 
([5]) for the DCE resonance with properly adjusted frequency shift, 77 = 2(cij — ^0)1 using the parameters rJo/'^ = 
1.4,^0/'^ = 2 • 10^^,£/ijj = 4 • 10^^. We show the number of created photons N (figured^) and the atom excitation 
probability Pg (figure [4|3, where the fast oscillations have been averaged out) versus the dimensionless time \80\ t for 
the initial ground atomic state \g) and the field number states |no) (no = 0, 1) and coherent states \g, a) (|a| = 1/2, 1). 
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For initial times, while (n) ric ^ 400, there is photon generation and amplification, as expected. However, the 
atom gets entangled with the field and grows as time goes on, suffering fast oscillations (not shown for clarity), as 
expected from equation ((35|) . 

In figures [3J; and HJl we show N and Pg, respectively, for larger times and different resonance shifts ^ = xJq, 
X — {0.95, 1, 1.05, 1.15, 1.25}, considering the initial state |.g,0). In the first place, we see that the photon generation 
is very sensitive to the resonance shift in the vicinity of ^ = (5, and near the 'resonance', a: « 1, there is a substantial 
photon creation from vacuum. However, the photon generation process ceases after several photons have been created 
and the photon number decreases, although the photon generation process restarts afterwards (data not shown). Pf. 
resembles the behavior of A^, but remains small for all the times. We also verified numerically that similar results are 
obtained for different atomic initial states, in agreement with equation ([55)1 . Therefore, for 77 = 2 (w ± (5) one obtains 
photon generation from vacuum, as well as photon number amplification, analogously to the dynamical Casimir 
effect 



VI. DISCUSSION OF RESULTS FOR PERIODIC MODULATION 



In previous sections we have considered only the first order resonances. In general, the K-ih order resonances occurs 
for an integer K when 



{K) _ j^-i 



(38) 



where ry^ stands for the AJC, JC and DCE resonances. In this case one recovers the previous results upon substitutions 



6 k and 



(39) 



as can be deduced from equations ^ and ([7]). For instance, to obtain 9k, one expands both sums on right-hand 
side of equation ([7]) and make RWA, keeping the terms which oscillate with the lowest frequency 1^. So 6k contains 
contributions due to the non-harmonic shape of the modulation (e.g. Ak e/ij), as well as due to the powers of {e/i]) 
[e.g. (Ai e/77)^ /i^!]. The effective dispersive shift Sk is obtained analogously, by first performing RWA in ([7]) and 
keeping terms which oscillate with frequencies higher than ^, and then calculating the effective Hamiltonian, as in 
equations (fMl) and So 5k contains the contribution of many terms besides 5q, among them the Bloch-Ziegert 

shift g'^/Aj^ [66] and powers of (s/rj). As an example, in figure we show N versus \0k\ got for the initial number 
state [17, 4) and the first three JC resonances, considering the resonance shift ^ = —IOSq. This example illustrates that 
higher order resonances are readily accessible. However, for the higher order resonances the effective dispersive shift 
6k should be carefully evaluated either numerically or experimentally, otherwise there is a risk of missing the exact 
resonance, since \9k\ becomes smaller and there is less freedom in committing small errors in the resonance shift ^. 

Our results can be directly transposed to the situation where n{t) — Q is constant and go = go (t) has a periodic 
time modulation, 



go {t) = go + eft- 

In this case, in the interaction picture the interaction Hamiltonian is 



Hr = 



{go + eco) - k {Ale 



ikr]t 



A.e-^'^^') 



k=l 



(40) 



(41) 



If one expands exp(iE!-|-) in equation ([7]) up to the first order in e/77, Hamiltonian ([5]) becomes equivalent to Hamil- 
tonian (|4T|) . so the results obtained above for ft {t) also hold for go {t) after making appropriate substitutions. The 
main difference between these two cases is that in the il {t) case the higher order K-th resonances occur both due 
to the powers of e/rj and non-zero coefficients Ak, while in the go {t) case the higher order resonances are due only 
to non-zero coefficients Ak- Therefore, for a harmonic modulation, ft = sinryf, in the fl{t) case one has higher 
order resonances, but in the go {t) case there are only the first order resonances. Finally, if the cavity frequency uj is 
modulated periodically, with constant fl and go, the AJC and JC resonances also occur, besides the well known DCE 
resonance [62, 63, 64, 65]. 

An experimental verification of this scheme seems possible in circuit QED architecture with superconducting qubits 
and coplanar waveguide resonators [38| , where one can adjust the system parameters in situ via electric and magnetic 
fields, as demonstrated in several experiments [ssl IstI l43j . Moreover, several schemes to read the cavity and the 
atomic states are currently available [sl, [3^ 113] (or under investigation ) . The main task would be finding a way 
to modulate periodically the atomic transition frequency with a stable modulation frequency i] 10 GHz, what is also 
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FIG. 5: The mean number of created photons N and atomic population inversion (sz) obtained numerically for the exact 
[equation (|3}] and effective [equation (|46p ] Hamiltonians for the resonance shift ^ — —Sq. The exact A*' and (az) show rapid 
oscillations, while the effective ones are smooth. 



within an experimental reach [37[. One could also use this scheme to couple M identical qubits (e.g. superconducting 
2-level atoms [13, El] or a cloud of polar molecules [s^, HH ) to the same cavity mode and modulate the frequency of 
AI atoms simultaneously, since in this case the effective coupling is increased to \/Mgo. 

One important point we did not analyze here is the dissipation and decoherence of both the artificial atom and 
the cavity due to the noisy solid state environment (48l . \6T\ . Recent experiments achieved experimental values 
{k/uj < 10^^,7/Li; < 10^'^,7p^/cLi < 10^^} [3§|, where k is the cavity decay rate, 7 is the atomic decay rate and 
7p;j is the atomic pure dephasing rate. To deal with dissipation in a qualitative manner, we compare the rates of 
the photon generation from vacuum for each resonance to the dissipation rates. We take the current experimental 
value of the coupling constant go/^ ~ 2 • 10~^ [s^ and assume e/u) ^ A^/lu ^ 10^^ to make the estimative. The 
photon creation rate for the first order DCE resonance is roughly \6d\ jio ~ 10^^ [equation ([35)) ]. and for the first 
order AJC resonance is roughly go \^\ ^ 10"'^ [equation (HH])]. Both these values are larger than the dissipation and 
decoherence rates (or of the same order of magnitude). Therefore, the photon creation due to modulation of 51 (t) or 
f/o(i) seems possible in the future. 



VII. A GLANCE AT NONPERIODIC MODULATION 



Let us now consider that the coupling parameter g (t) has a bichromatic external modulation 

9 (i) = .90 + 2£sin[(A+ - 2^)^] + 2£_ sin[(A_ - 2^_)t], (42) 

where go is the bare coupling parameter, \e\ , ]e_] <C 50 are small modulation amplitudes and ^, ^_ are the 'resonance 
shifts'. Here A± is given by equation (jTU]) with cq = (i.e., A± = ± u) and we assume A_ > 0. The system 
Hamiltonian is given by equation ([3]) with time-dependent g (t) and constant lu and Q. In the dispersive regime, we 
obtain the effective Hamiltonian 

HeffC^6{h+ 1/2) a, + {iee^'^*a^a+ + ie-e^'^-*aa+ + h.c). (43) 

Performing the unitary time-dependent transformation 

C/3 = exp [(e + C-) <^z/2 + (e - C-) n] t} (44) 
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we get the time-independent Hamiltonian combining the 'rotating' and 'anti-rotating' terms with adjustable couphngs 

Hnp ^ C/3 HeffUs - zt/g 

C + +(5(l + 2n) , N / + N , ^ 

= (7z + — n + [iea^a+ + ie-aa+ + h.c.) . (45) 

We take equal resonance shifts in order to cancel the second term in (|45|) . £,=£,_, so for £_ = (e = 0) we obtain 
the effective Anti-Jaynes-Cummings (Jaynes-Cummings) Hamiltonian. A more interesting regime occurs when both 
e and e_ are different from zero. For e = £_ we obtain 



^ -cr^+ie{a + a'^){cr+-cr^), (46) 

while for e = — £„ we get pS)) with the last term replaced by ie (a — a^) (ct+ + (t_). The Hamiltonian cannot 
be integrated exactly, so below we solve it numerically and compare the results to the ones obtained using the initial 
Rabi Hamiltonian ([3]). 

To obtain photon creation from the initial vacuum state \g, 0), we have to make the exponent multiplying in (|46p 
as small as possible by adjusting the resonance shift ^. In figure [5] we show the mean number of created photons N 
and the population inversion (sz) as function of dimensionless time et for e = S- and ^ = — i5o, calculated numerically 
for the initial Hamiltonian ([3]) and the effective Hamiltonian (|46p using the parameters ^l/cu = 1.3, go/ui = 5 • 10~^ 
and e/u! = 5 • 10"'^. One can see that the effective Hamiltonian describes well the exact dynamics, and there is 
a quasi-periodic generation of a few photons, together with the atomic excitation. Qualitatively, this occurs because 
for initial times the Hamiltonian becomes 

H+ ~ -~V2exay, (47) 

and the photons are generated from vacuum quadratically in time. However, after 1 photons have been created, 
the exponent multiplying az in (j46|) becomes large, so the photon creation process goes out of resonance and it is 
interrupted. 




FIG. 6: The mean number of created photons A'^ obtained numerically for the exact [equation ((3|] and effective [equation (|46|l ] 
Hamiltonians for the resonance shift £, = £, (t) [equation (148 The effective A'^ shows slow oscillations, while the exact one does 
not. 

To create a large number of photons from vacuum, we can adopt an 'active' approach, by which the resonance shift 
f is adjusted adiabatically as function of time, so the effective Hamiltonian (|45p is valid. In this way, we continuously 
adjust the resonance shift to the number of photons in the cavity, provided we choose an appropriate functional form 
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of ^. In figure[6]we show the mean number of created photons A'' from the initial state \g, 0) for the Rabi Hamiltonian 
([3]) and the effective Hamiltonian using the parameters Vt/uj = 1.4, go/uj — 5 • 10~^, e/u) = 5 • 10^^ and setting 

i^l{t) = -{5^~ 3.g„VA+) (1/2 + et) . (48) 

One can see that a significant amount of photons can be generated from vacuum, provided the resonance shift is 
adjusted adiabatically, and the effective Hamiltonian describes well the exact dynamics. 

Finally, notice that these results also hold if, instead of the coupling parameter one modulates the atomic 
transition frequency fJ, or both g and simultaneously, one with frequency ~ A-|_, and another with ~ A_. 

VIII. CONCLUSIONS 

In conclusion, we studied the nonstationary solid state cavity QED system in which the atomic transition frequency 
or the atom-cavity coupling have a small periodic time modulation, prescribed externally. We showed that in the 
dispersive regime and under the resonant modulations, the Rabi Hamiltonian (which describes the system dynamics) 
can be significantly simplified, resulting in three different regimes that can be described by the Anti-Jaynes-Cummings, 
Jaynes-Cummings or the dynamical Casimir effect Hamiltonians. Moreover, in the resonant atom-cavity regime, 
entangled states with two photons can be created from the vacuum state |(7,0) under the corresponding resonance, 
analogously to the dynamical Casimir effect in a vibrating cavity containing a resonant two-level atom. 

We also showed that the photon generation from vacuum occurs for a small nonperiodic time modulation of the 
atom-cavity coupling, the atomic transition frequency, or both in circuit QED. We deduced an effective Hamiltonian 
for the bichromatic modulation in the dispersive regime, and demonstrated that it describes well the exact dynam- 
ics. Numerical simulations confirmed that several photons can be generated from vacuum provided the modulation 
frequencies are fine tuned adiabatically. 

This study illustrates the importance of the anti-rotating terms in the Rabi Hamiltonian, ignored in the Jaynes- 
Cummings model - here this term is responsible for photon generation from vacuum and field amplification. Moreover, 
one could engineer effective interactions in nonstationary circuit QED by means of modulating the system parameters. 
As applications, this scheme can be used to verify photon creation from vacuum in nonstationary circuit QED due 
to an analog of the dynamical Casimir effect using a single atom, as well as off-resonant population transfer between 
the states {\g,n), |e,rt± 1)} and generation of entangled states with several photons. 
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